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Abstract 

We present an improved measurement of the branching fraction for the electroweak penguin 
process B Xsi^i~, where i is an electron or a muon and Xg is a hadronic system containing 
an s-quark. The measurement is based on a sample of 152 x 10® T(4S') ^ BB events collected 
with the Belle detector at the KEKB energy asymmetric e'’“e“ collider. The Xg hadronic system 
is reconstructed from one or and up to four pions, where at most one pion can be neutral. 
Averaging over both lepton flavors, the inclusive branching fraction is measured to he B{B 
Xgi+e-) = (4.11 ±0.83(stat)lg:g(syst)) x 10“® for M^+^- > 0.2 GeV/c^. 

TAGS numbers: 13.20.He, 12.15.Ji, 14.65.Fy, 14.40.Nd 


*on leave from Nova Gorica Polytechnic, Nova Gorica 
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I. INTRODUCTION 


In the Standard Model (SM), the rare decay B —»• proceeds through a 6 —>• 

transition, which is forbidden at tree level. Such a flavor-changing neutral current (FCNC) 
process can occur at higher order via electroweak penguin and W^W~ box diagrams. The 
b —> transition therefore allows deeper insight into the effective Hamiltonian that 

describes FCNC processes and is sensitive to the effects of non-SM physics that may enter 
these loops; see, for example. Refs 0,1^. 

Recent SM calculations of the inclusive B —> branching fractions for Mi+i- > 0.2 

GeV/c2 predict B(B X^e+e") = (4.2 ± 0.7) x lO-^ 0, 0, and B{B X,£+£-) = 
(4.6 ±0.8) X 10“® 0. Both the Belle and BaBar collaborations have observed exclusive B 

and decays 13, ra, ItI and have measured the rate for inclusive B —>• Xsi^i~ 

decay Hi. 

In this analysis, we study the inclusive B —> Xsi^i~ process by semi-inclusively recon¬ 
structing the hnal state from a pair of electrons or muons and a hadronic system consisting 
of one or Xg and up to four pions, where at most one pion can be neutral. This semi- 
inclusive reconstruction approach US allows approximately 53% of the full inclusive rate 
to be reconstructed. If the contribution of the modes containing a is taken to be equal 
to that containing a Xg, the missing states that remain unaccounted for represent ~30% 
of the total rate. We require the hadronic mass for the selected hnal states to be less than 
2.0 GeV/c^ to reduce combinatorial background. We correct for the missing modes and the 
effect of the hadronic mass requirement to extract the inclusive B —>■ Xsi~^i~ decay rate 
for Mi+i- > 0.2 GeV/c^. This measurement updates and supersedes our previous result 0, 
which was based on a sample of 65 x 10® BB pairs. 


II. THE BELLE DETECTOR AND DATA SAMPLE 

We use a data sample collected on the T(4S') resonance with the Belle detector at the 
KERB energy asymmetric e’''e“ collider (3.5 GeV on 8 GeV) HH. This sample comprises 
152 X 10® B meson pairs, corresponding to an integrated luminosity of 140 fb“^. A detailed 
description of the Belle detector can be found elsewhere ^2. A three-layer silicon vertex 
detector (SVD) and a 50-layer central drift chamber (GDG) are used for tracking and iden- 
tihcation of charged particles. An array of aerogel threshold Cerenkov counters (AGG) and 
time-of-hight scintillation counters (TOF) are also used for charged particle identihcation. 
An electromagnetic calorimeter comprised of Tl-doped Gsl crystals (EGL) measures the en¬ 
ergy of electromagnetic particles and is also used for electron identihcation. These detectors 
are located inside a superconducting solenoid coil that provides a 1.5 T magnetic held. An 
iron hux-return located outside of the coil is instrumented with resistive plate counters to 
identify muons (KLM). 

Particle identihcation for e^, X^, Xg, and vr® is important for this analysis. Elec¬ 

tron identihcation is based on the ratio of the cluster energy to the track momentum [E/p), 
the specihc energy loss {dE/dx) in the GDG, the position and shower shape of the cluster in 
the EGL and the response from the AGG. Muon identihcation is based on the hit positions 
and the depth of penetration into the EGL and KLM. Electrons and muons are required to 
have laboratory-frame momenta greater than 0.4 GeV/c and 0.8 GeV/c, respectively. To 
select good muon candidates, we apply a kaon veto. We hnd an electron (muon) selection 
efficiency of 94.3% (93.2%) with a 0.1% (0.9%) pion to electron (muon) mis-identihcation 
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probability. Bremsstrahlung photons from electrons are recovered by combining each elec¬ 
tron with photons within a small angular region around the electron direction. Charged 
kaon candidates are selected by using information from the ACC, TOF and CDC. The kaon 
selection efficiency is 90% with a pion to kaon mis-identification probability of 6 %. After 
selecting the electron, muon and charged kaon candidate tracks, the remaining tracks are 
assumed to be charged pious. ATg candidates are reconstructed from pairs of oppositely- 
charged tracks with \M^+^- — M^o| < 15 MeV/c^ (~ Sctm q). We impose additional iCg 

S -KTg 

selection criteria based on the distance and the direction of the Kg vertex, and on the impact 
parameters of the daughter tracks. We require all charged tracks, except those used in the 
Kg reconstruction, to have impact parameters with respect to the nominal interaction point 
of less than 1.0 cm in the radial direction and 5.0 cm along the beam direction. Neutral pi¬ 
ous are required to have a laboratory-frame energy greater than 400 MeV, photon daughter 
energies greater than 50 MeV, and a 77 invariant mass satisfying \M^^ — M^o \ < 10 MeV/c^ 


III. ANALYSIS OVERVIEW 


We reconstruct inclusive B —>■ decays with a semi-inclusive reconstruction tech¬ 

nique from a pair of electrons or muons and one of 18 different hadronic states. Here the 
hadronic system consists of one or Kg and up to four pions (at most one pion can be 
neutral). Compared to a fully inclusive approach, this method has the advantage of having 
strong kinematical discrimination against background by using the beam-energy constrained 
mass Mbc = \J— p\ and the energy difference KE = Eb — Abeam, where Abeam is the 
beam energy and Eb {pb) is the reconstructed B meson energy (3-momentum). All quan¬ 
tities are evaluated in the e^e“ center-of-mass system (CM). 

Further background suppression to reduce the large combinatorial backgrounds is neces¬ 
sary. The main contribution to the combinatorial background comes from pairs of semilep- 
tonic 6 —>• c decays in A A events. In these events, A —> Xsi^i~ candidates are reconstructed 
with the decay products from both AA mesons. This background has a significant amount 
of missing energy due to the neutrinos from the semileptonic decays. Another contribu¬ 
tion to the combinatorial background comes from continuum events, which are efficiently 
suppressed with event-shape variables. 

There are two background sources that can peak in Mbc and AA. The first comes 
from the charmonium decays A — J/xjjX and A — 'ip{2S)X with J/'ijj{'ip{2S)) —>■ Ki~. 
This charmonium background is efficiently removed with requirements on the dilepton mass 
Mi+i-. The resulting veto sample provides a large control sample of decays with a signature 
identical to that of the signal. The second comes from A — K^{Kg)n7r [n > 1) decays with 
misidentification of two charged pions as leptons. We estimate these peaking background 
contaminations, then subtract them from the signal yield. 

For the A —>■ Xsi^i~ event simulation, we use EVTGEN Q as the event generator, JET- 
SET Q to hadronize the system consisting of a strange quark and a spectator quark, and 
GEANT for the detector simulation. In the event generation, A —>■ Xsi~^i~ events are 
produced with a combination of exclusive and inclusive models. In the hadronic mass region 
Mxs < IT GeV/c^, exclusive A —^ decays are generated according to Refs |lL 

where the relevant form factors are computed using light-cone QGD sum rules. In the region 
Mxg > 1.1 GeV/c^, event generation is based on a non-resonant model following Refs (lLll7j| 
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and the Fermi motion model of Ref. Hi. 


IV. EVENT SELECTION 

Events are reqnired to have a well determined primary vertex, be tagged as mnlti-hadron 
events, and contain two electrons (or mnons) having laboratory-frame momenta greater than 
0.4 GeV/c (or 0.8 GeV/c for mnons). Dilepton candidates consist of or pairs. 

The two leptons are reqnired to originate from a common vertex and satisfy the reqnirement 
I A^l < 0.015 cm. Here, Az is the distance between the two leptons along the beam direction; 
the z-coordinate of each lepton is determined at the point of closest approach to the beam 
axis. 

Gharmoninm backgronnds are rednced by removing B candidates with a dilepton mass in 
the ranges —0.40 GeV/c^ < Mee(^) — Mj/^ < 0.15 GeV/c^, —0.25 GeV/c^ < — Mj/^ < 

0.10 GeV/c2, -0.25 GeV/c^ < M^( 2 S) < 0.10 GeV/c^, and -0.15 GeV/c^ < 

M^( 2 S) < 0.10 GeV/c^. If one of the electrons from a J/'0 or decay erroneonsly picks 

up a random photon in the bremsstrahlung-recovery process, the dilepton mass can increase 
sufficiently to evade the above requirements. Therefore the charmonium veto is applied to the 
dilepton mass before and after bremsstrahlung recovery. Using the simulation, we estimate 
the remaining peaking charmonium background to be 1.20 ± 0.28 events and 1.33 ± 0.21 
events for the e^e“ and channels, respectively. 

The requirement Mi+i- > 0.2 GeV/c^ is applied. For the e+e“ channel, this removes 
potential peaking background from the conversion of photons from radiative B —^ 
decays and from 7r° Dalitz decays. 

Starting from a pair, B —>■ candidates are formed by adding either 

a or a Kg and up to four pious, but no more than one 7r°. In this manner, 
eighteen different hadronic topologies are reconstructed: K^, K^tt^, K^tt'^tt^, 

K^TT^TT^, iF^TT^TT^TT^, iF^TT^TT^TT^, iF^TT^TT^Tr^Tr^, iF^Vr^TT^TT^Tr^, Kg, KgTT^ , KgTl^ , 
KgTT^TT^, KgTT^TT^, iFgTT^Tr^Tr^, iFgTT^Tr^Tr^, iFgVr^Tr'^Tr^Tr^, and iFgVr^Tr^Tr'^Tr^. 

After forming the B —>• Xs£~^£~ candidates, we apply additional requirements to suppress 
background. The largest background sources are random combinations from continuum qq 
{q = u, d, s, c) production or from semileptonic B decays. We reject the qq background by 
a requirement on a Fisher discriminant [l^ (Afw) based on a modihed set of Fox-Wolfram 
moments Q that characterize the event topology. In the semileptonic B decay background, 
both B mesons decay into leptons or two leptons are produced from the b ^ c ^ s,d 
decay chain. This B decay background is rejected by the requirement on another Fisher 
discriminant variable (Amiss) constructed from the total visible energy Avis = Ei and the 

missing mass Mmiss = \J (2Abeam — S EiY — \ I]Pip, where (p), Ei) are the reconstructed 
four-momenta, in the GM frame, of all tracks (assumed to be pions) and all photons in the 
event. 

We further reduce the background by requiring —0.10 GeV < AA < 0.05 GeV for 
the dielectron channel (—0.05 GeV < A A < 0.05 GeV for the dimuon channel), and 
xlt^/NDF < 10. Here, Xvtx is the of fh® B vertex constructed from the charged daugh¬ 
ter particles, excluding the Kg daughters. We also reject candidates with an X^ invariant 
mass greater than 2 GeV/A. This condition removes a large fraction of the combinatorial 
background while retaining 99% of the signal. 

At this stage, there is an average of 1.6 A candidates per event in the signal simulation. 
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In order to select the most signal-like B candidate, we construct a likelihood ratio based on 
the following six discriminant variables: Ai?, = -Eroe — -E'beam, Xvtx; cos 9 b, ^fw, 

and j Fmis s. The energy of the rest of the event, -Eroe is calculated by summing the energies 
of all charged tracks and neutral calorimeter clusters not included in the B candidate, and 
cos 6 b is the cosine of the B flight direction with respect to the e~ beam direction in the 
CM frame. The variables AE, AE^^^, and are effective at rejecting BB background, 
especially for events with two semileptonic decays, which have large missing energy. For 
continuum suppression, the event-shape variables Epw and cos 9 b are useful. The variable 
is effective for rejecting the random combinatorial background in the high multiplicity 
modes. 

The probability density functions (PDFs) of these six variables are determined from MC 
events, separately for signal and background, with the exception of signal PDFs for AE and 
Xvtx- determine the latter two from the charmonium data sample, which satishes all 
selection criteria except for the charmonium veto (charmonium-veto sample), because we 
observe some discrepancy between the signal MC and the charmonium-veto data sample. 
To determine the PDF from data, we subtract the background using the shapes obtained 
from the Mbc sidebands. The number of background events in the Mbc signal region is 
normalized by htting a Gaussian (signal) and an ARGUS function (background) to the 
Mbc distribution. For AE, we use separate PDFs for the dielectron and dimuon channels. 

We then calculate the likelihoods £sig and Drg which are products of the PDFs of the 
above six discriminants for the signal and the background, respectively. In each event, only 
the B candidate with the largest value of the likelihood ratio IZ = £sig/(£sig + 'Cbg) is 
retained. According to the MG simulation, this choice implies that 84% of the selected B — > 
candidates in R —> events have all daughter particles correctly assigned. 

To check the distributions of the likelihood ratio 71, we compare data and MG, separately 
for signal (Fig. ^) and background (Fig. ^) events. For the signal, we use charmonium- 
veto candidates from data and B —>■ Xsi^i~ signal MG simulation. In the 7Z distribution 
for charmonium-veto events, we subtract the background using shapes obtained from the 
Mbc sidebands and normalization from the number of background events in the Mbc signal 
region, in the same way as for the signal PDF determination. For the background we use 
B —>■ candidates that are reconstructed using the nominal selection criteria but 

requiring that the two leptons have different flavor, in data and in a background MG sample 
that contains 6 — »• c decays and qq events. We observe good agreement in both cases. 

The hnal suppression of the combinatorial background is achieved with a requirement 
on the likelihood ratio 7Z, optimized to maximize the statistical signihcance of the signal 
S/\/S + B. Here S and B denote the numbers of the signal and background events, re¬ 
spectively. This MG optimization is performed separately in the regions Mx^ < 1.1 GeV/c^ 
and 1.1 GeV/c^ < Mxs < 2.0 GeV/c^, resulting in the requirements 7Z > 0.3 and 7Z > 0.9 
respectively. 

After applying all selection criteria, 155 B —> Age+e" and 112 B As/i+p“ candidates 
remain in the Mbc signal region, dehned as 5.27 GeV/c^ < Mbc < 5.29 GeV/c^. According 
to the simulation, the background remaining at this stage of the analysis consists mostly 
of BB events (80% and 72% of the total background in the electron and muon channels, 
respectively). In MG signal events, the probability that all daughter particles from a selected 
B —> As£+£“ candidate are correctly assigned is 91%. 
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FIG. 1: Likelihood ratio for (a) signal and (b) background events. Points represent data and solid 
histograms represent MC. Samples used are described in the text. 


V. MAXIMUM LIKELIHOOD FIT 

We perform an extended unbinned maximum likelihood fit to the Mbc distribution in the 
region Mbc > 5.2 GeV/c^ to extract the signal yield as well as the shape and yield of the 
combinatorial background. The £t is performed separately for the dielectron and dimuon 
channels. The likelihood function C is expressed as: 

p-N N 

■c=irrr ni(«'.ig+iVpe.k)pr* ++ Navf + iVc.„bP“”'’i 

i=i 


where N = Wig + Weak + We + Wf + Womb IS the total number of candidate events. Wig, 
Weak, Wc) Wf, and Womb represent the yields of the signal, peaking background, combina¬ 
torial background from peaking background, cross-feed background (i.e. events containing 
a true B decay, but a selected B candidate with incomplete or wrong decay 

products) and combinatorial background, respectively. As explained below, we use the same 
PDF, Pf®, for the signal and peaking background components. 

The signal PDF Pf® is described by a Gaussian for the channel as well as for 

the channel, since the bremsstrahlung recovery and selection procedure for the 

channel lead to a negligible radiative tail in the Mbc distribution. The Gaussian shape 
parameters are determined from fits of the sum of a Gaussian and an ARGUS function 
to the charmonium-veto data sample. The hts result in a mean Mbc of rrisig = 5279.31 ± 
0.05 MeV/c^ and a Mbc resolution of cTsig = 2.62 ± 0.04 MeV/c^ for the e+e“ channel, and 
rusig = 5279.03 ± 0.04 MeV/c^ and dgig = 2.53 ± 0.04 MeV/c^ for the channel. In 

the simulation, the Gaussian fit results for the Mbc distributions of correctly reconstructed 
signal events are in agreement with the shape parameters extracted from the hts to the 
charmonium-veto sample. The signal yield Wig is a free parameter in the likelihood £t. 

We consider two separate contributions to the peaking background: the charmonium 











(a)Xs charmonium : e+e- (b) Xs charmonium ; |li+|li- 



(c)Xs h+h-: e+e- (d) Xs h+h-: |li+|li- 



FIG. 2: Mbc distributions for MC charmonium events for the (a) dielectron and (b) dimuon 
channels, and for B —> Xsh'^h~ candidates in data for the (c) dielectron and (d) dimuon channels. 
The normalization corresponds to the expected charmonium and hadronic peaking background in 
140 fb“^ of data. 


peaking background, already described, and the hadronic peaking background. The latter 
arises from the B —> (n > 0) —> Xs7r~^Tr~ decays where two hadrons are misidentihed 

as leptons. The normalization and shape of this background are determined directly from 
the data by repeating the selection without the lepton identihcation requirements, and 
htting the Mbc distribution with the sum of a Gaussian and an ARGUS function. We 
obtain rupeak = 5279.16 ± 0.04 MeV/c^ and Upeak = 2.60 ± 0.02 MeV/c^, consistent with 
the signal shape of the charmonium peaking background. Hence we use a common PDF 
for these two peaking backgrounds, which is the same as the one determined for the signal. 
Taking into account the momentum and angular dependence of the tt —> £ misidentihcation 
probabilities (on average 0.08% for electrons and 0.92% for muons) the hadronic peaking 
background is estimated to be W-peak = 0.030 ± 0.001 events for the dielectron sample and 
A^h-peak = 1-78 ± 0.05 eveuts for the dimuon sample. The charmonium peaking background 
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TABLE I: Results of the fit to the data: number of signal candidates in the Mbc signal box, fitted 
signal yield Ngig, estimated peaking backgrounds Apeak (fixed in the fit), and significance including 
systematics. 


Mode 

Candidates 

^Sig 

■^peak 

Significance 

A, e+e- 

155 

31.8 ± 10.2 

1.24 ±0.28 

3.2 

As 

112 

36.3 ± 9.3 

3.11 ±0.22 

4.4 

Combined 

267 

68.4 ± 13.8 

4.35 ± 0.36 

5.4 


is estimated from the simulation to be Ac_peak = 1-20 ± 0.28 and iVc_peak = 1-33 ± 0.21, 
respectively. In the likelihood hts, the peaking background normalization is fixed to Apeak = 
Ac_peak + Ah_peak- FiguresEfa) and (b) show the Mbc distributions for the MC charmonium 
events for the dielectron and dimuon channels, and (c) and (d) show the Mbc distributions 
for B —> Xsh^h~ candidates in data for the dielectron and dimuon channels, respectively. 

The other background PDFs, Vf, and are given by an ARGUS shape. They 

describe the combinatorial contribution from peaking background events, from cross-feed 
events, and from continuum and BB events, respectively. The ARGUS cutoff is determined 
by the beam energy in the T(4S') rest frame, Abeam = 5.290 GeV. The values of the AR¬ 
GUS shape parameter for each background component are determined from the peaking 
background Mbc distribution shown in Fig. El from incorrectly reconstructed signal 

MG events {Vf), and from the fit to the B —> data selected using the nominal 

selection criteria but requiring that the two leptons have different flavor 
these three ARGUS shape parameters. We also fix the number of peaking background events 
Ape and the number of cross-feed events Acf, from the peaking background Mbc distribution 
shown in Fig. Eland from incorrectly reconstructed signal MG events, re^ectively. Here the 
signal MG generations are based on the branching fractions in the Ref. [l| . The yield Acomb 
is taken as a free parameter in the likelihood fit. 


VI. RESULTS 

Using the fit parametrization described above, we fit the Mbc distributions in data sep¬ 
arately for the B Xse^e~ and B —>■ The fit results are displayed in Fig. Eland 

summarized in Table HI The significance is 5 = 2 Amaij, where £max represents 

the maximum likelihood for the £t and denotes the maximum likelihood when the 

signal yield is constrained to be 0. To include the systematic uncertainty of the signal yield 
estimation into the signihcance, one of the htting parameters, which are signal Gaussian 
parameters (mean and width), the yield of the peaking BG, and the ARGUS shape param¬ 
eters for the background PDFs is varied within Tier in the above signihcance calculation. 
We then take the minimum significance value. The B —> Xs^'^^~ signal yield presented in 
Table n is obtained by a £t to the combined electron and muon samples. Figure EKd) shows 
the Mbc distribution for B —>• Xge^jjX candidates. Applying the ARGUS £t to the Mbc 
distribution, there is no evidence for a peaking background as expected. 

Figures EKa) and (b) show the distributions of the hadronic mass Mx^ and = M|b£- 
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(a)Xse+e- (b)Xsn+H- 




FIG. 3: Mbc distributions of (a) B Xse~^e~, (b) B —> (c) B —> Xsi'^i~ {I = e, ^), and 

(d) B —> XsC^ijX candidates in data. The solid lines represent the result of the fits, and dashed 
lines represent the sum of all background components under the signal peaks, respectively. 


for the electron and muon channels combined, obtained by repeating the likelihood £t in 
bins of Mxs and Figure EJ^a) indicates that the observed signal includes contributions 
from hnal states across a range of hadronic masses, including hadronic systems with a mass 
above that of the K*{892). 

The branching fraction B for the signal is calculated as 


B 




( 1 ) 


where = (152.0 ± 0.7) x 10® is the number of BB pairs produced in 140 fb ^ and e is 
the signal efficiency. 


VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are of two different types: those that affect the extraction of the 
number of signal events and those that affect the calculation of the branching fraction. The 
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FIG. 4: Distributions for the signal yield as a function of (a) hadronic mass Mx^ and (b) 

for B —> signal. The points with error bars are the data (electron and muon channels 

combined) while the histogram is the MC signal normalized to the data. 


systematic uncertainties are summarized in Table UTl 

Uncertainties affecting the extraction of the signal yield are evaluated by varying the sig¬ 
nal Gaussian parameters (mean and width) and the background shape parameter within Tin 
of the measured values from the charmonium-veto data sample (signal) and the B —> 
data sample (background). We estimate the uncertainties in the peaking-background shape 
and cross-feed events by comparing the signal yields obtained with and without the corre¬ 
sponding PDFs in the unbinned maximum likelihood £t to Mbc- 

Uncertainties affecting the signal efficiency originate from the detector modeling, from the 
simulation of signal decays, and from the estimate of the number of B mesons in the sample. 
By far the largest component is that due to the simulation of signal decays, discussed in 
detail below. 

The detector modeling uncertainty is composed of the following uncertainties determined 
from the data: the uncertainty in the tracking efficiency of 1.0% per track; the uncertainty 
in the charged-particle identihcation efficiency of 2.05% per electron, 2.95% per muon, 1.0% 
per kaon and 0.8% per pion; and the uncertainty in the reconstruction efficiency of 4.5% per 
Kg and 3.3% per vr®. The efficiency of the likelihood ratio requirement, which suppresses 
combinatorial background, is checked with the charmonium-veto sample and the level of 
discrepancy with the simulation is taken as the corresponding uncertainty. 

The dominant source of uncertainty arises from modeling the signal decays. Parameters 
of the Fermi motion model are varied in accordance with measurements of hadronic moments 
in semileptonic B decays and the photon spectrum in inclusive B ^ Xg j decays [ 2 ^ . 
The fractions of exclusive B —>■ Ki~^£~ and B —>■ K*Ki~ decays are varied according to 
theoretical uncertainties P|. The transition point in Mx^ between pure K*Ki~ and non¬ 
resonant Xg Kl~ final states is varied by ±0.1 GeV/c^. 

The non-resonant Monte Garlo event generator relies on JETSET to fragment and 
hadronize the system consisting of a final state s quark and a spectator quark from the 
B meson. Since the signal efficiencies depend strongly on the particle content of the £- 
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TABLE II: Relative systematic uncertainties (in %) on the B Xge~^e and B Xg^'^ 
branching fraction measurements. The different contributions are summed in quadrature. 


Source 

X,e+e- X, 


Signal shape 

±1.4 

±0.5 

BG shape 

±7.8 

±4.7 

Peaking background statistics 

±0.9 

±0.6 

Peaking background PID error 

< 0.1 

±0.5 

Peaking background shape 

±4.3 

±2.1 

Cross-feed events 

±4.1 

±2.2 

Signal yield total 

±9.9 

±5.7 

Tracking efficiency 

±3.5 

±3.5 

Lepton identification efficiency 

±4.1 

±5.9 

Kaon identification efficiency 

±0.8 

±0.8 

TT^ identification efficiency 

±0.6 

±0.5 

Kg efficiency 

±0.7 

±0.8 

TT® efficiency 

±0.3 

±0.3 

7Z requirement efficiency 

±5.4 

±4.5 

Fermi motion model 

- 1 - 6.5 

- 2.4 

-re.i 

- 2.3 

B(B K£+£-) 

±9.9 

±10.5 

B(B K*£+£-) 

±7.0 

±7.8 

K*-Xs transition 

±4.5 

±4.7 

Hadronization 

±8.5 

±8.2 

Missing modes 

±4.5 

±4.4 

Monte Carlo statistics 

±1.6 

±1.5 

Efficiency total 

- 1 - 19.0 

- 18.1 

- 119.7 

- 18.9 

BB counting 

±0.5 

±0.5 

Total 

- 1 - 21.5 

- 20.6 

- 120.5 

- 19.7 


nal state, uncertainties in the number of charged and neutral pions and in the number of 
charged and neutral kaons translate into a signihcant uncertainty in the signal efficiency (for 
Mx, >1.1 GeV/c2). 

The ratio of the generator yield for decay modes containing a Kg to that for modes 
containing a charged kaon is varied according to 0.50 ±0.11, to allow for isospin violation 
in the decay chain. The ratio of the generator yield for decay modes containing one 
meson to that for modes containing none is varied according to 1.00 ± 0.22. Uncertainties 
in these two ratios are set by the level of discrepancy between B — J/ipX data and MC 
B Xgi~^i~ samples. 

The 18 modes selected in this analysis only capture about 53% of the full inclusive rate. 
If the contribution of the modes containing a is taken to be equal to that containing a 
Kg, the missing states that remain unaccounted for represent ~30% of the total rate. To 
account for the uncertainty in such remaining missing states, we estimate the uncertainty in 
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TABLE III: Summary of results: signal yield (Ngig), statistical significance, total signal efficiency 
e (including the fraction of Xg states considered in this analysis) and branching fraction (B). The 
first and second errors quoted on Ngig and B are statistical and systematic, respectively. The 
first error on e corresponds to uncertainties in detector modeling, BB counting, and Monte Carlo 
statistics, and the second error on e to the uncertainties in the signal model. 


Mode 

Nsig 

Significance 

6(%) 

B (xlO-6) 

Xg e~^e~ 
Xg 

A, £+£- 

31.8 ± 10.2 ±3.1 

36.3 ± 9.3 ±2.1 

68.4 ± 13.8 ±5.0 

3.2 

4.4 

5.4 

2.59 ±0.20 toil 
2.89 ±0.24:toil 
2.74 ±0.22lg-j| 

4.04 ± 1.30 

4.13 ± 1.05 Igjf 
4.11 ±0.83lg;|f 


TABLE IV: B ^ Xgt^l branching fractions (B) for different bins in Mx^ and The 

first and second errors are statistical and systematic, respectively. 


Mxs in GeV/c^ 

B (xlO-^) 

in (GeV/c)^ 

B (xlO-^) 

[0.4, 0.6] 

[0.6,0.8] 

[0.8,1.0] 

[1.0,1.4] 

[1.4, 2.0] 

3.75 ± 0.96l[jJ^ 
0.36 ± 0.88l[j;[}^ 
6.65 ± 2.25^0^7 
10.50 ± 6.90l|[]^ 
46.59 ± 23.37tllil 

[0.04,1.0] 

[1.0, 6.0] 

[6.0,14.4] 
[14.4,25.0] 

11.34 ± 4.83l|^? 
14.93 ± 5.04l|^} 
7.32 ± 6.141};®^ 
4.18±1.17l°;6i 


the fraction of modes with too many pions or kaons (two extra kaons may be produced via 
ss popping), as well as the contributions of modes with photons that do not originate from 
7r° decays but rather from rj, r/', etc. For hnal states with Mx^ > 1-1 GeV/c^, we vary these 
fractions by ±5% per 7r°, ±20% for r/, ±30% for > 5, and ±50% for rj' and others. 

Including systematic uncertainties, the measured branching fractions for M^+i- > 
0.2 GeV/c^ are 

B{B X,e+e-) = (4.04 ± 1.30x 10“®, (2) 

B{B Xgjj+fi-) = (4.13 ± 1.05+°!®) X 10“®, (3) 

B{B^Xgi^r) = (4.11 ±0.83l|]|f) X 10“®, (4) 

where the hrst error is statistical and the second error is systematic. The combined B —> 

Xg£~^i~ branching fraction is the weighted average of the branching fractions for the electron 
and muon channels, assuming the individual branching fractions to be equal for M^+i- > 
0.2GeV/cT Table Hm summarizes the results of the analysis and lists both the statistical 
and systematic errors on the signal yields, the signal efficiencies and the branching fractions. 

The branching fractions for each Mx^ and bin are also measured, and summarized in 
Table HVl Figures EKa) and (b) show the distributions of the differential branching fractions 
as a function of hadronic mass Mx^ and respectively, for electron and muon 

channels combined. 
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FIG. 5: Differential branching fraction as a function of (a) hadronic mass and (b) 
for B —> signal. The points and error bars represent the data (electron and muon channels 

combined) while the histograms represent the MC signal normalized to the data statistics. The 
outer (inner) error bars represent the total (statistical) errors. 


VIII. SUMMARY 

Using a sample of 152 x 10® T(4S') —> BB events, we measure the branching frac¬ 
tion for the rare decay B —^ where i = e or fi and Xg is a hadronic system 

reconstructed semi-inclusively in 18 different hadronic states (with up to four pions). For 
Mi+i- > 0.2 GeV/c^, we observe a signal of 68.4 ± 13.8{stat) ± 5.0{syst) events and obtain 
a branching fraction (averaged over lepton species) of 

B{B Xgtr) = (4.11 ± 0.83(stat) ^a8?(syst)) x 10“®, 
with a statistical signihcance of 5.4 a. 

This result is consistent with the recent prediction by Ali et al. [H 0, 0| , our previous 
inclusive B —> Xgt^t~ measurement [^, and that of the BaBar collaboration j^. 
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